Astronomy & Astrophysics manuscript no. nh ©ESO 2015 

February 27, 2015 


Long term variability of Cygnus X-1: 

VII. Orbital variability of the focussed wind in Cyg X-1 / HDE 226868 system 

V. Grinber^, M.A. Leuteneggei®!], N. Heli?0, K. Pottschmid01l, M. Bocll^ J.A. GarcftP, M. Hanke!^, 
M.A. Nowal!® J.O. Sundqvis£0, R.H.D. TownsencP^, and J. WilmP 


in 

H 

O 

(N 

X) 

D 

IJ-H 


' Massachusetts Institute of Technology, Kavli Institute for Astrophysics and Space Research, Cambridge, MA 02139, USA 
^ CRESST, University of Maryland Baltimore County, 1000 Hilltop Circle, Baltimore, MD 21250, USA 
^ NASA Goddard Space Flight Center, Astrophysics Science Division, Greenbelt, MD 20771, USA 
Dr. Karl-Remeis-Stemwarte and ECAP, FAU Erlangen-Nurnberg, Sternwartstr. 7, 96049 Bamberg, Germany 
^ Lawrence Livermore National Laboratory, 7000 East Ave., Livermore, CA 94550, USA 
® Max-Planck-Institut fur Radioastronomie, Auf dem Htigel 69, 53121 Bonn, Germany 
’ Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA 02138, USA 
* University of Delaware, Bartol Research Institute, Newark, DE 19716, USA 

® Institut fur Astronomie und Astrophysik der Universitat Munchen, Scheinerstr. 1, 81679 Munchen, Germany 
Department of Astronomy, University of Wisconsin-Madison, 2535 Sterling Hall, 475 North Charter Street, Madison, WI 53706, 
USA 


Received: — / Accepted: — 

(N 



Oh 

6 

a 


ABSTRACT 

Binary systems with an accreting compact object are a unique chance to investigate the strong, clumpy, line-driven winds of early type 
supergiants by using the compact object’s X-rays to probe the wind structure. We analyze the two-component wind of HDE 226868, 
the 09.7Iab giant companion of the black hole Cyg X-1 using 4.77 Ms of RXTE observations of the system taken over the course of 
16 years. Absorption changes strongly over the 5.6 d binary orbit, but also shows a large scatter at a given orbital phase, especially 
at superior conjunction. The orbital variability is most prominent when the black hole is in the hard X-ray state. Our data are poorer 
for the intermediate and soft state, but show signs for orbital variability of the absorption column in the intermediate state. We 
quantitatively compare the data in the hard state to a toy model of a focussed Castor-Abbott-Klein-wind: as it does not incorporate 
clumping, the model does not describe the observations well. A qualitative comparison to a simplified simulation of clumpy winds 
with spherical clumps shows good agreement in the distribution of the equivalent hydrogen column density for models with a porosity 
length on the order of the stellar radius at inferior conjunction; we conjecture that the deviations between data and model at superior 
conjunction could be either due to lack of a focussed wind component in the model or a more complicated clump structure. 
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1. Introduction 

Early type supergiants show strong line driven winds (CAK 
mechanism; [Cast or, A bbott, & KleTn| |1975t |Morton| |1967[ 
|Lucy & Solomon|[lj70^nd typical mass loss rates of a few 


10 ”Mo/yr (e.g.. 


Puls et al. 


2006| l. These high velocity winds 
(e.g.,|Mi 


that can reach terminal velocities Vcx, > 20 00 km s _ 

jres et al.|2012|l are perturbed and clumpy (|Owocki et al.|1988 


Feldmeier et al. 1997 


Sundqvist & Owocki 


Send offprint requests to: V. Grinberg, 
e-mail: grinberg@space.mit.edu 


Puls et aL]|20061 lOskinova et al.||2012 
2013| l with over 90% of the wind mass 
concentrated in less than ~10% of the wind volume (|Sako et al.| 
[1999) [Rahoui et al.)|20TT) for Vela X-1 and Cyg X-1, respec¬ 
tively). Binary systems consisting of an O/B type giant and a 
compact object offer us the unique chance to investigate these 
winds by using the X-ray source as a test probe. 

Cygnus X-1 is in such a high mass X-ray binary with 
the 09.7Iab supergiant HDE 226868 ( Walborn]|1973|l. L ocated 
~1.86kpc away ( |Reid et ar]|2011| see also Xiang et ar]|201 l| l, 
it is one of the brightest and best-observed black hole bina¬ 
ries. The orbital period is P = 5.599829 + 0.000016 d with 


To = 52872.288 + 0.009 ( |Gies et al.||2008] l. The source also 
shows a superorbital period that has been reported as ~300d 
(|Priedhorsky et al.|1983 1, then as ~150d ( Brocksopp et al.|1999| 


[Benlloch et al.|2004[[Ibragimov et al.|2007|l, and lately as having 


changed back to ~300 d (jZdziarski et al.| 26 lT i. 


The intrinsic variability of Cyg X-1 on timescales of hours 


(|Bock et al.|201 1 

\ to years (jPottschmidt et al. 2003| Wilms et al.| 

2006| Shaposhni 

cov & Titarchuk 

2006 Axelsson et al. 2006) 

Grinberg et al.|2013|l can be classified in terms of X-ray spectral 


states that also show distinct timing characteristics on timescales 
below 1 s and are similar for most black hole binaries (for a re¬ 
view see, e.g., |Belloni|2010| ). In particular, the hard state spec¬ 
trum above ~2 keV is dominated by a power law component with 
photon index T ~ 1.7, while in the soft state the thermal emis¬ 
sion from the accretion disk is prominent and the contribution 
from a steeper power law component is low. 

jOrosz et al.j ( |201 l| l determine the inclination of the binary 
system to i - 27?1 ± 0?8, the black hole mass to Mbh = 
(14.8+1.0)Mo, and the companion mass to M* = (19.2+1.9)Mo, 
the latter value being questioned by Ziolkowski ( 2014| l, who esti¬ 
mate a range of 25-35 Mq with a most likely value of 27 Mq. As 
HDE 226868 fills >90% of its Roche lobe ( |Gies et al.|200^ and 
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has a strong (a few 10 ® Mq yr ', 

Gies et al. 200311 wind, the 

black hole accretes via a focussed wind ([Friend & Castor 

1982 

Gies & Bolton[[1986a|bl Sowers et al.|19981 Miller et al. 

2005 

Hanke et al.[[2009 Hell et al. 20131 Miskovicova et al. 

2015 

and others). Values for for HDE 226868 vary in the litera- 


ture: 


Davis & Hartmann (|1983 1 obtain (2300 + 400) km s 


state. 


the h ard state. butIVrtilek et al. (|2008ll 1420 km s ' in the soft 

' jGies et~ar| ( |2008 


Herrero et al. 


(1995 I adopt 2100kms 


as 


obtain 1200kms ^ during a soft state, but list 1600kms 
more likely. 

The inclination of the system implies that our line of sight 
probes different regions of the wind at different orbital phases. 
Orbital variability of the absorption column density in the hard 
state has been reported by Feng & Cui| ( |2002 1, [Ibragimov et al.| 
( 2005| l, and Boroson & Vrtilek ( 2010| ). Dips in the (soft) X-ray 
lightcurves of Cyg X-1 are interpreted as signatures of clumps 
in the stellar wind. They occur preferably at superior conjunc¬ 
tion (0orb ~ 0), i.e., when the O-star companion is between the 
observer and the black hole (e.g., M & Clark] 1974l|Mason et al. 


1974| Pa-signault et al.||1976 jPravdo et al.||1980 Remillard & 

Caniz^es|1984[|Kitamoto et al.|1984[|Balucinska-Church et al. 


2000 Poutanen et al.|2008 Boroson & Vrtilek|2010 l. 

The 16 years of mainly bi-weekly Cyg X-1 observations 
(1996-2011) with RXTE’s Proportional Counter Array (PCA; 
Jahoda et al.j 2006|l and High Energy X-ray Timing Experi¬ 


ment (HEXTE; Rothschild et al.|[l998| l allowed us to study the 
long-term spectral and X-ray timing evolution of the source 
in unprecedented detail in the previous paper of the series 


([Pottschmidt et al. [[20031 

Gleissner et al.[[2004b|al [Wilms et al.[ 

[2006 

Grinberg et al. 

201 

3 

2014|l. In this paper, we use these 


data to analyze the orbital variability of absorption due to the 
focussed wind. The exceptional exposure reveals subtle effects 
such as the orbital variability in intermediate states and rare 
events with very high absorption. 

In Sect.j^ we introduce the data used and the spectral models 
applied to describe them. In Sect. we discuss the observed 
orbital variability of absorption in hard, intermediate, and soft 
states. In Sect. we compare our measurements with previous 
results and with different wind models, namely with a toy model 
for a CAK-wind and with simulations of clumpy winds in single 
stars. We summarize our results in Sect. [5] 


2. Data and spectral analysis 


[Ibragimov et al.|f2005| l used 42 Ginga-OSSE and RXTE-OSSE 
observations of Cyg X-1 taken mainly in the hard state to show 
orbital variability of absorption by demonstrating an increased 
equivalent hydrogen column density, Ah, around 0orb ~ 0. In 
[Wilms et al.[ ( 2006| l, we analyzed 202 observations from 1999- 
2004 using RXTE-PCA and, where available, -HEXTE spectra 
with models similar to the one presented below and with a time 
resolution of individual RXTE observations that were typically 
of several satellite orbits in length, i.e., had exposure times of 
5-10 ks. [Boroson & Vrtilek[ ( [2010[ ) classified 102 of these ob¬ 
servations as hard states and used them to demonstrate the or¬ 
bital variability of absorption in the hard state, but did not have 
enough data taken during softer states. At the same time, higher 
sensitivity snapshots at different (porb with Chandra ([Miller et al.[ 


[20051 [Hanke et al.|2009 

Hanke|20111 Miskovicova et al.|20151l 

and Suzaku ([Nowak et a 

i.|201 lllMiller et al. 2012[|Yamada et al.[ 


able local densities, but the orbital coverage was significantly 
worse than the RXTE measurements. 
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exposure time [ks] 


Fig. 1. Histogram of the PCA exposure times of the individual satellite 
orbit-wise observations analyzed in this work. The mean exposure is 
1.74 ks, the median 1.66ks. 



Fig. 2. Number of orbitwise PCA observations with a given soft pho¬ 
ton index F] that do (red) or do not (blue) require a multitemperature 
disk component. Additionally highlighted are observations that do not 
require a disk and fall into the long hard state of MJD 53900-55375 
(black). Dashed vertical lines indicate thresholds for hard (Fi < 2.0), 
intermediate (2.0 < Fi < 2.5) and soft (2.5 < Fi) states according to 
[Grinberg et al.[j2()13) . 


In this paper we therefore extend the earlier RXTE anal¬ 
yses, using the full mission-long database of RXTE observa¬ 
tions of Cyg X-1, and using spectra taken with the exposure of 
one orbit of the RXTE satellite around Earth, or around 1.6 ks 
(Fig. a- We thus probe average spectral parameters on shorter 
timescales than earlier RXTE studies over a variety of spec¬ 
tral shapes (Eig. |^. Our total sample consists of 2741 PCU 2 
and HEXTE-spectra ( [Grinberg et al.[[^13| l. The total exposure 
time of our observations is 4.77 Ms which corresponds to al¬ 
most 10 fully covered orbital periods of 5.6 d. The exposure is 
evenly distributed between the orbital phases. Advances in PCA 
calibration ( [Shaposhnikov et ah 20121 and a careful consider¬ 
ation of energy-dependent systematic uncertainties in different 
gain epochs of the PCA-instrument (Hanke[201 1[ Grinberg et al. 


2013[l strengthen our analysis. Eor a detailed discussion of data 


treatment, see [Grinberg et al.[ ( [2013[ l. All analyses in this pa 
per were performed with ISIS 1.6.2 ([Houck & Denicola|[20()0 
[Houck|2002[ [Noble & Nowak|2008[ ). 

In performing our analysis we used the same approach as 
that employed by us previously ( Grinberg et al.[[20T3 l, i.e., we 
model the data empirically with a broken power law with soft 
photon index Ei, hard photon index r 2 , and a spectral break at 
~10keV ( Wilms et al. [[2006 1. This broken power law is modi¬ 
fied by a high-energy cutoff, by an Ee Ko'-line modelled with 
a Gaussian at ~6.4keV, and by absorption. A multitemperature 
disk (diskbb, [Mitsuda et al.[|T984| [Makishima et al.[|1986[ l is 
added to the model where it improves the by more than 5% 
except in seven cases where the X-ray timing behavior and the 
correlation between Fi and F 2 strongly prefer the model without 
a disk ([Grinberg et al.|20T4|). Figure]^ shows that the model of- 
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fers a very good description of the spectrum. As in |Grinberg et al.| 
( 2013| l, we use Fi-based state definitions derived from spectral 
and timing properties of Cyg X-1; the source is in hard state if 
Fi < 2.0, in the intermediate state if 2.0 < Fi < 2.5 and in the 
soft state if 2.5 < F]. The distribution of observations with Fi is 
shown in Fig.[^ Out of 1822 hard state observations, 1515 do not 
require a disk, but 384 out of 413 intermediate state observations 
and all 506 soft state observations require a disk component. 

When studying variations of Ah in an astronomical source, 
the measured column towards the source consists of two main 
components; absorption in the interstellar medium (ISM) that 
is constant on observable timescales and a variable absorption 
local to the X-ray source. For Galactic sources. X-ray data do 
not allow us to distinguish the two absorption components. In the 
case of Cyg X-1, we therefore take the constant ISM contribution 
into account by setting a min imum value for Ah to Ah.ism = 
4.8 X 10^' cm“^, as derived by Xiang et al. ( 201 1| ) from studies 
of the dust scattering halo of the source. We describe abs orption 
with the tbnew model, an improved versiorpjof tbabs (Wilms 


et al.|2000|l, using wilm abundances of Wilms et al. (|2060|l aiid 


vern cross sections of |Verner et al.| ( |i996| ). Note that optical 
depth at a frequency v depends on the total cross-section cry and 
the hydrogen column density Ah as Ty — cry ■ Ah. The observed 
spectrum, /obs, and the source spectrum, 4^, are then related by 


^obs.i 


= L, 


For the individual contribution to cTy from 


different phases of the absorbing medium (e.g., gas or dust) and 
the different constituents (e.g., oxygen) see [Wilms et al.|p000| l 
or jGatuzz et al.| ( |2()15| l. Gener ally in the X-ray ra nge, the optical 
depth is dominated by metals ( Wilms et al.|2000 \ and cry oc 
i.e., the photoelectric absorption acts more strongly on the low 
energy photons. 

We caution that the exact choice of absorption model, abun¬ 
dances, and cross-sections can lead to differences in Ah up to 
20%-30% between different approaches ( Wilms et al.|2000| , al¬ 
though they do not affect oyerall trends obseryed. We also cau¬ 
tion that the absorption modeling is based on the absorption of 
neutral material, while in reality the material will be moderately 
photoionized by the strong X-ray source embedded in it. Espe¬ 
cially when using proportional counters such as the PCUs, ap¬ 
proximating the X-ray absorption by ionized material with ab¬ 
sorption in neutral material is an appropriate assumption, since 
most absorption is due to K-shell absorption, which is only 
mildly dependent on ionization stage. Howeyer, our modelling 
will not pick up effects from fully ionized material, as it is mainly 
transparent to X-rays and RXTE’s instruments did not haye the 
resolution to detect absorption lines. 

A second issue relates to the uncertainty of the Ah measure¬ 
ments. Eor models such as the ones used here there is a well 
known degeneracy between the power law slope, i.e., Fi, and 
Ah (e.g., jSuchy et al.||200^ . By calculating two-dimensional 
confidence contours we confirm that this degeneracy is not re¬ 
sponsible for the yariation of Ah discussed here (for example 
contours see |Hanke|201 Ij l. Because of the limited energy range 
of the PCA, in addition to the correlation with the power law 
continuum there is also a degeneracy between parameters of the 
accretion disk and Ah. The effect of a stronger absorption can be 
counterweighted by a stronger disk component and yice yersa. 
An-yalues from best fit models with and without a disk therefore 
haye to be considered separately because of the different system¬ 
atic effects and we do so throughout this paper. Unless discussed 
otherwise, we giye all uncertainties at the 90% leyel. 


* http: //pulsar. sternwarte. uni - er langen. de/wi 1ms/ 
research/tbabs/ 



Fig. 3. Histogram of values for observations that do (red) and do 
not (blue) require a multitemperature disk component. 
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Fig. 4. Example spectra with 1.72 < Fi < 1.73, but different absorption 
columns. PCA data and best fit models to PCA and HEXTE data are 
shown. The ratio plots show residuals of the best fit in color and resid¬ 
uals of an evaluation of the best fit model with Ah = Ah.ism (i-e., no 
absorption local to the system) in black. All observations are from the 
hard state of 2006-2010 (Sect. 

3. Orbital variability of absorption 

3.1. Hard state 

We now turn to studying the variability of Ah in the RXTE 
data. As an example. Fig. shows spectra of different obser¬ 
vations which have the same typical hard state continuum shape 
(Fi = 1.72 ... 1.73) but different absorption columns. The resid¬ 
uals shown in the figure for various orbital phases illustrate the 
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Fig. 5. Absorption column vs. unabsorbed energy flux during the hard 
state of 2006-2010. Red circles mark observations shown in Fig.[^ 


strong effect caused by the orbital variation compared to the in¬ 
fluence of the interstellar medium alone. Absorption affects the 
whole spectrum up to energies of lOkeV. 

When looking at the whole state-resolved sample, we deter¬ 
mine the average and median values of Ah and its standard de¬ 
viation, cr{AH, 22 ) := cr{AH/(10^^ cm“^)), for phase bins of width 
A0orb = 0.1 starting at (^oi-b = 0. In the hard state, the average col¬ 
umn values and the standard deviation, i.e., the variability of the 
wind between individual observations, increase around superior 
conjunction (Fig.[^ left). 

Individual observations, especially for 0orb ~ 0, show Ah 
values much larger than the average for a given phase (Fig. 
left, see also Fig. for example spectra). We associate these 
observations with the prolonged deep absorption events (often 
dubbed “dips”) known from the literature (Sect. and care¬ 
fully check our data to find out whether these observations have 
systematic differences to other observations in the same spectral 
state. We find no such differences; in particular, these “dips” do 
not have larger uncertainties in Ah and do not show clear trends 
in their exposure time, unabsorbed 3-50 keV flu)0(see also 
Fig. or the correlation between Fi and F 2 known from|Wilmk 


et al. 


P006| l. The dips cluster around superior conjunction, but 
individual strong absorption measurements occur as early as at 
0orb ~ 0.8 and as late as at (^orb ~ 0.3 (Fig.[^ left). 

The RXTE spectra do not allow a meaningful fit of partial 
covering models and so we cannot constrain possible partial cov¬ 
ering either in physical space or in time (i.e., a clump or a group 
of clumps, possibly of different optical depth, passing through 
the line of sight for a part of the total exposure time of one spec¬ 
trum). Thus, the An-values derived for individual observations 
are lower limits on the sum of the absorption columns of individ¬ 
ual clumps with the ISM absorption over a given observation’s 
exposure time (i.e., on average 1.6-1.7 ks, Fig.[2l. 

Within the hard state, Cyg X-1 can display a range of spec¬ 
tral and timing behaviors ([Pottschmidt et al.|2003[ [Wilms et al. 


2006 1 Grinberg et al. 2014| l. At the same time, we expect the 
behavior of the absorption to depend on the broadband spectral 
shape of the X-ray source because of possible ionization of the 
wind. For a detailed analysis of absorption, we therefore need a 
hard state with stable spectral and timing characteristics, such as 
the long hard state of 2006-2010 (MJD 53900-55375, jGrinbergj 


et al. 


riod in Fig. 
a disk (Fig. 


2bl3||20T4|. We show results for this long hard state pe- 


1\ it includes 383 observations that do not require 
2l and that have a total exposure time of 660 ks or 


-1.4 orbital periods of the system. The variability of absorption 


^ The X-ray flux of Cyg X-1 can vary by factor of 4-5 at the same 
spectral shape 1 Wilms et al.|2006 and references therein). 
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can be clearly seen in the (cumulative) probability distributions 
of Ah values at different orbital phases (Fig.[^ left). The values 
are low and the distribution narrow at 0orb ~ 0.5. At (^oib ~ 0, the 
distribution is broad with a higher average (average values and 
their standard deviations are shown in Fig.[^ right). We will use 
this set for comparisons with a toy model for the focussed wind 
and a clumpy wind model in Sect.|4.1.^and|4.13] respectively. 


3.2. Intermediate state and soft state 


We now address observations of Cyg X-1 in the intermediate and 
the soft state, i.e., when the X-ray spectrum is much softer than 
in the hard state. Most data in the intermediate state require a 
disk component that leads to larger uncertainties in the Ah deter¬ 
mination and particularly to a systematically increased average 
value independent of orbital phase as well as a larger cr{AH, 22 )- 
An increase of average Ah at <porh ~ 0 is, however, still visible 
(Fig.|^ middle). 

Our data do not show orbital variability of absorption in the 
soft state (Fig.[^ right). The average values are comparable with 
the intermediate state, but in several phase bins we find median 
values of 4.8x10^' cm“^, i.e., the data are consistent with no neu¬ 
tral absorption intrinsic to the source. Because of the increased 
disk contribution, observations in the soft state are even more af¬ 
fected by systematics than those in the intermediate state, espe¬ 
cially given that the disk in Cyg X-1 has a temperature < 0.4 keV 
(=s 4.6 X 10® K) even in the soft state ( Wilms et al.|2006l and our 
fits), i.e., it peaks below the lower limit of our data so that the 
disk parameters show relatively large uncertainties. 


4. Discussion 

4.1. Hard State 

4.1.1. Discussion and comparison to earlier results 

Our results show that despite PCA’s energy coverage only above 
~3keV and comparatively low spectral resolution when com¬ 
pared to CCD or even grating instruments, PCA data can be used 
to assess trends in the behavior of the absorption. The changes in 
Ah strongly influence the shape of the observed spectra (Fig.|^. 
The overall orbital variability of Ah values found here is consis¬ 
tent with earlier results (Feng & Cui 2002 [Ibragimov et al.|2005j 


[Wilms et al.[2006[[Boroson & Vrtilek[20ib| l. To our knowledge 

we are the first to quantify the variance of the equivalent hydro¬ 
gen column, cr{AH, 22 ) (Fig-|^ lower panels). 

The largest Ah values seen during the hard state (Figs. 
and 1^ are on the order of lO^^cm^^. ||Kitamoto et al. (1984 


using Tenma), [Balucinska-Church et ah ( 1997[ using ASCA), 


and |Hanke et al. ( 2008[ using Chandra) reported individual deep 


absorption events with Ah ^ lO^^cm Periods of increased 
absorption that lasted longer than a few ks have been observed 
with various instruments, especially close to superior conjunc¬ 
tion. For example, Feng & Cui ( 2002| l discuss a long dip of at 
least ~3ks that shows complex substructure. Prolonged periods 
of increased absorption are also clearly seen in the long uninter¬ 
rupted observations possible with Chandra ( [Hanke et al.|2008} . 
Given these previous measurements and the size of our sample, 
the detection of several strong absorption events as discussed in 
Sect. |3.1| is plausible. 

[Balucinska-Church et al. (2000i have used mainly RXTE- 
All Sky Monitor (ASM; Levine et a .|1996| ) data to show orbital 
variability in ASM-based hardness ratios and in dip (defined by 
a threshold ASM hardness ratio) occurrence, both peaking at su- 
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orbital phase 0orbit 



orbital phase 0orbit 



b -0.5 0 0.5 1 1.5 


orbital phase 0orbit 


Fig. 6. Variation of Nu with orbital phase for observations where the source was in the hard (left), intermediate (middle), and soft (right) state. The 
upper panels show the orbit dependence of Nu measurements for individual satellite orbit-wise spectra (black datapoints with gray uncertainties), 
as well as the average (red) and median (white) values for orbital phase bins with a width of A(f>orb = 0.1. The lower panels show the standard 
deviation cr{A^H 22 l •= cr|^H/(10^^ cm“^)) of the measurements. The panels specifically show 1515 hard state (Ti < 2.0) observations that do 
not require a disk component, 384 intermediate state (2.0 < Ti < 2.5) observations that require a disk, and 506 soft state (2.5 < Ti) observations 
(all soft state observations require a disk component). 



Fig. 7. Observations from the hard state of 2006-2010 (MJD 53900-55375, [Grinberg et al.|2013t that do not require a disk. Left: normalized 
histograms (blue) and normalized cumulative histograms (red) for Nh above the ISM value for individual orbital phase bins; n denotes the total 
number of observations in a given orbital phase bin. Note that the phase bin ^ = 0.7-0.8 contains one measurement in excess of 10^*’ cm“^. Right: 
Orbital dependence of Nu- The best fit focussed wind model introduced in Sect.|4.1.^is shown in green. 


peri or conjunction consistent with the behavior of the absorption 
column density that we observe. They observe a peak in the dip 
occurence at (^orb = 0.95, but our sample is not large enough 
to track this possible assymetry of the A^h distribution around 
0orb = 0. [Poutanen et al.j ( |2008| l further show that the orbital 
modulation of ASM countrates in the hard state depends on the 
superorbital phase: the number of ASM-defined dips shows a su¬ 
perorbital periodicity when assuming a ~150d superorbital pe¬ 
riod. The smaller number of our measurements when compared 
to ASM data combined with the apparent stochastic nature and 
large variance of the absorption variations do not allow such an 
analysis here. 

Despite these shortcomings, the exceptional coverage of the 
pointed observation allows us to compare the observational re¬ 
sults to theoretical models. In particular we quantitatively com¬ 
pare a simple focussed CAK wind to our data from the stable 
hard state period of 2006-2010 (~250 orbits of the binary sys¬ 
tem) in Sect. [4.1.2| and, for the first time for a high mass black 
hole binary, qualitatively match the observed data from the same 
stable hard state period with predictions from a clumpy wind 
model for O-star winds in Sect. 14.1.31 


4.1.2. Focussed wind model in the hard state 


We first compare the orbital variation of Ah from MJD 53900- 
55370 with a toy model for the focussed wind as introduced by 
[Gies & Boltonj (|1986b|l. T his model is based on the previous 
work of Friend & Castorj (1982l and depends on the angle 6 
between the wind direction and the donor-black hole axis. We 
emphasize that this is only a toy model that does not take into 
account details of both the physics of a real wind and the ob¬ 
servations used. The model is, however, useful to assess general 
trends and to point out the directions in which better models have 
to be developed. 

For 0 ^ 20°, the wind is the radiatively driven (CAK-)wind 
of a single star ( [Castor et al.|1975] l. For 0 < 20°, the wind varies 
smoothly with 0 with terminal velocity and density peaking at 
0 = 0. In particular, we employ the model with the fill-out-factor 
of 0.98 with the correspondent values for terminal velocity and 
wind density listed in Table 1 of [Gies & Bolton ( 1986a|l a nd use 
masses, inclination, and orbital period as given in Sect.[]J Figure 
l^shows the resulting distribution of wind density. 
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X 


Fig. 8. Stellar wind density calculated following the focussed wind 
model of |Gies & BoltonH1986a| ). The x- and z- axes are shown in units 
of orbital separation. The position of the black hole is marked with a 
black cross, the position of the center of mass with a green cross. Lines 
of sight towards the black hole at orbital phases 0 and 0.5 are indicated 
by dashed lines, the orbital plane by a dotted line. Note that there are 
no discotinuities in the wind; if such appear above, they are unfortunate 
tricks of vision. 


We obtain the total (neutral and ionized) column density 
by integrating along the line of sight towards the black hole 
at a given orbital phase: these values are on the order of 3- 
5x lO^^cm^^, i.e., much higher than the measured average values 
and more reminiscent of dips (cf. Fig. left) But this total col¬ 
umn density is not the A^h measured in our observations: the fully 
ionized part of the wind is transparent to X-rays and thus only 
some weakly to moderately ionized part of the wind will con¬ 
tribute to the measured absorption. We can hence assume that 
the absorption in the binary is restricted to the region of the wind 
wi th log f < log frtian with the ionization parameter f defined af¬ 
ter 


Tarter et al. 


(1969 1 as ^ = Ll(nn) with n being the absorbing 


particle number density, r the distance from the ionizing source, 
and L the luminosity above the hydrogen Lyman edge (ty pical 

We 


Miskoyicoya et al 


dgejtyp 


for Cyg X-1: L - 10 ergs 
then consider only the contribution of those parts of the wind 
that satisfy this condition. 

We take into account the ISM column density A^h.ism and 
fit the parameter log^max of this model to the data. We ob¬ 
tain log^max ~ 2.7, howeyer with = 9.1 for 382 degrees 
of freedom. The phase yariation of A^h that corresponds to our 
best fit yalue of log ^niax is shown in Fig. right panel, as a 
green curye. The failure to describe the data well is expected 
since the toy model does not include clumping that we expect 
to majorly contribute to the spread of the A^h measurements at a 
giyen orbital phase. To conduct a simple test for the influence of 
clumping on our results, we remoye the 13 measurements with 
A^h >5x 10^^ cm“^ and obtain a bestfit;i^J^j = 4.5 for 369 degrees 
of freedom with a decreased bestfit yalue (log^max ~ 2.67 as 
compared to preyiously log^max ~ 2.71). 

The toy model also neglects additional smaller effects such 
as a potential photoionization wake that would introduce asym¬ 
metries ( |Blondin||1994| l. Uncertainties in determination of the 
ISM absorption, A^h.ism, towards the source may lead to fur¬ 
ther problems ( Xiang et al.||201l) . Howeyer, it is interesting to 
note that when we account for the fact that our measurements 
include the clumping, our result is in rough agreement with re- 
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suits (log^ < 2.6) found by Miskoyicoya et al. ( 2015 1, who infer 
A^h yalues ranging from 6.8 x 10^^ cm“^ to 1.3 xTO^ cm“^ from 
modelling the Ne edge in the non-dip parts of hard state Chandra 
obseryations at 5 different orbital phases. 


4.1.3. Clumpy wind model for the hard state 


We next consider whether the obseryed yariability in absorbing 
column could be due to wind structures, and how the character¬ 
istics of those structures compare to what is known about radia- 
tiyely driyen winds in other contexts, such as clumps ([Abbott 


eTaflfTWTl |Hillier & MiUer][TW8l |Fullerton et ar]|2006t [Puls' 


et al.|2006l or discrete absorption components (DACs; Prinja & 


Howarth||1986|l. To this end, we use a numerical code to sim¬ 


ulate the column density between the obseryer and a point-like 
X-ray source, as the source orbits around a companion whose 
wind is composed of discrete clumps. For simplicity, the code 
does not account for the focussed wind properties of the sys¬ 
tem; moreoyer, the clumps are assumed to be spherical, rather 
than the pancake morphologies considered by other authors (e.g., 
Feldmeier et al.|2003[ Oskinoya et al.|200j|. At a giyen orbital 


phase, the column-density eyaluation is performed by casting a 
ray from the source to the obseryer, determining which ray seg¬ 
ments lie inside clumps, and summing up the contributions from 
these segments (as the product of the segment length and the 
clump density). 

The distribution of clumps throughout the wind follows a 
similar approach to that described in Appendix A of|Sundqyist 


|et al.| ( |2OT^ . Clumps are relased at the stellar surface r = /?, in 
a random direction and with an initial radial size /», and then ad- 
yect outward with the wind according to a canonical ‘/I’ yelocity 
law v(r) - Voo(l-Rtlr'f. As a clump adyects its radial size grows 
as l{r) = {rlRt)l,, but its mass remains fixed at m - M/N where 
M is the oyerall wind mass loss rate, and N is the clump re¬ 
lease rate. The latter is most conyeniently parameterized in terms 
of the number of clumps per wind flow time N - iVfflow, with 
fflow = R*lvca- To relate these parameters to the porosity length h 
(essentially, the mean free path between clumps; see |Owocki &| 
|Cohen|2006| l, we note that h(r) - haai>(r)IVoc, where the terminal 
porisity length is 


hoc = 3- 


R, 


( 1 ) 


To apply the code to Cyg X-1, we assume p - 1.6, = 

I^Rq, Vcc - 2400kms“*, and M = 2 x lO^^Moyr^', and 
adopt the orbital parameters giyen in Section 1. We ran a grid 
of simulations with N - 100,316,1000,3162 and 10000, and 
hoaIR* = 0.1,0.316,1,3.162, and 10 to explore the leyel of yari¬ 
ability expected on the time scales we probe here. For each grid 
point we calculated 1000 model orbits, with a phase resolution 
A0oib = 0.001. Since typical RXTE obseryations of Cyg X-1 had 
an exposure of about 0.005 orbits, and since each obseryation is 
not subdiyided into time slices, it is necessary to ayerage the 
model to a timescale corresponding to the obseryations. Thus, 
we rebinned the models by a factor of fiye by taking a linear ay¬ 
erage of the model column density. A linear ayerage oyer a phase 
bin does not strictly correctly account for the non-linear addition 
of partial coyering in time, but is adequate giyen the crude na¬ 
ture of our model. Likewise, rebinning by a factor of fiye is not 
an exact match to the actual obseryations, but it is an adequate 
approximation. 

The statistical properties of each phase bin of 0.005 orbits do 
not change rapidly as a function of orbital phase, so these statisti¬ 
cal properties were further ayeraged into phase bins of 0.1 orbits. 
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Fig. 9. Predictions of A^h variability assuming a clumpy wind model with the total number of clumps per flow time N = 1000 and the terminal 
porosity length hoo = R, (see Sect. |4.1.3 for explanation of parameters). Left: model histograms for the distribution of equivalent absorption column 
density A^h. Right: average values (circles) and standard deviations (‘error bars’ on the average values) in a given orbital phase for the theoretical 
model (red) and values measured in the hard state of MJD 53900-55375 (black). Average measured values are shown as ‘average - A^h.ism’ to 
account for the ISM absorption. 


corresponding to the grouping used to evaluate the observations. 
We show this statistical analysis of the column for N = 1000 
and hoc - R* in Fig. the mean and standard deviation of this 
model are compared to those of the data in the right panel. We 
have chosen a model with a similar level of variability as the ob¬ 
servations near phase 0.5, where our line of sight to the compact 
object mainly probes the part of the wind that is likely least per¬ 
turbed by the peculiarities of the focussed part of the wind. For 
other model realizations see Figs.|A.l]and|A.2| 


The choice of porosity length hoc has a strong effect in deter¬ 
mining the level of variability in the model (Fig. |A.l| l. The choice 
of N has a weaker effect, with larger values of N (and thus small 
values of /») showing somewhat reduced variability (Fig. |A.2| i. 
This is a result of single observations averaging over the transit 
of one or more clumps across the line of sight, thus washing out 
the variability expected on short timescales. We have explored 
the effect of different choices of averaging timescale (or equiva¬ 
lently observation duration) in our simulations, and as expected, 
we find that increasing the averaging timescale tends to reduce 
the variability (see also Sect. 3.1 for a discussion of exposure 
times on the measured Ah). We also find that there is a maximum 
useful time resolution above which increasing the sampling rate 
does not increase the level of variability observed. This max¬ 
imum rate corresponds roughly to the time scale for a typical 
clump to transit our line of sight to the compact object which is 
on the order of k/voc (when the wind velocity dominates over 
orbital velocity as is the case here). One obvious conclusion of 
these results is that the clump size scale could be constrained 
by obtaining spectra of high statistical quality (e.g., with Suzaku 
or Chandra satellites) and exploring the level of variability on 
shorter time scales. 


The models with hoc - R* have the best agreement with the 
data near tpoib = 0.5, with the exception that for N - 10000 the 
model with hoc - 3.162/?* has the best agreement (Fig. |A.2] i. 
Based on previous theoretical and observational investigations 
of single O stars, it is reasonable that hoc should be comparable 
to or s maller than R, (|Dessart & 0wocki|2003||0wocki & Cohen| 
|2006||Naze et al.|20i3|l. Th is is also consistent with the results 

of Leutenegger et al. ( |2013| l, who found that the spectrum of the 
single O supergiant ( Pup could be fit with models having poros¬ 
ity lengths hoc < R*. To test for the influence of terminal velocity 
on our results, we re-ran the simulations with Voc - 1600km s“^. 
Although the resulting mean values for Ah are higher (due to the 
increased wind density), our results are otherwise similar to the 


original iioc - 2400 km s ' grid, and thus our conclusions do not 
change. 

Comparing Fig. and the model shows a Gaussian dis¬ 
tribution of absorbing columns, while the data show a non- 
Gaussian tail near phase 0. It seems likely that this non-Gaussian 
tail reflects some unusual wind structure related to the focussed 
wind. However, since the impact parameter of our line of sight 
with respect to the optical companion is smallest near this orbital 
phase, the tail might also be a signature of a wind with pancake¬ 
shaped clumps, which would be effectively more porous when 
viewed edge on ( [Feldmeier et al.|2003[|0skinova et al.|2012| l. A 
more thorough exploration of these effects may give insight into 
their importance but it is beyond the scope of this paper. 


4.2. Intermediate state and soft state 


Due to the lack of suitable observations in the previously ana¬ 
lyzed samples, the variability of absorption in intermediate and 
soft states has not been previously analyzed using RXTE spec¬ 
tra, but only using ASM data. Boroson & Vrtilek ( 2010| l ana¬ 
lyzed orbital variability of absorption in five soft state periods. 
[Boroson & Vrtilek| see no variability in their soft states 2 and 5 
that would be mostly classihed as intermediate states using the 
more rigorous state definitions of jGrinberg et al. ([2013 . Their 
soft state 3 shows orbital variability. Boroson & Vrtilek (|2010|l 


also point out that the intrinsic variability of the X-ray source in 
the softer periods hampers the search for periodicities in ASM 
data. 

Based on the pointed RXTE data alone, we see no orbital 
variability in the soft state, while there are signs for an increase 
of the main absorption toward (^orb ~ 0 in the intermediate state. 
These results qualitatively agree with the behavior seen in anal¬ 
yses of Chandra high resolution spectra: as the source softens, 
dipping becomes less pronounced ( jMiskovicova et sT. 2015| l. 
However, dips have still been observed in the soft state with 
Suzaku ( [Yamada et al. 2013| but note the strong changes in hard¬ 
ness ratio in their data that may be indicative of spectral evo¬ 
lution during the observation). Hanke| (20111 shows lightcurves 
and softness ratios of all Chandra observations up to and includ¬ 
ing ObsID 13219 taken in 2011 February; some weak and short 
dipping can be seen in the soft state data. 

The interpretation of these results in the context of the previ¬ 
ous Suzaku and Chandra results is that the soft X-rays emitted 
from the vicinity of the black hole in the soft state are strong 
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enough to fully ionize the wind and even some or all of the opti¬ 
cally thick clumps. The wind thus becomes (mainly) transparent 
to X-rays and no orbital variation of absorption column can be 
seen. 


5. Summary and conclusions 


We have analyzed almost 5 Ms of RXTE observations of Cyg 
X-1 taken over a span of 16 years in different spectral states. 
We find clear orbital variability of absorption in the hard state 
of the X-ray source with an overall increase and additional deep 
absorption events (“dips” with A^h ~ lO^^cm^^) around supe¬ 
rior conjunction, in agreement with previous work with different 
instruments. Because of the presence of an additional multitem¬ 
perature disk component in the spectra, the absorption column 
cannot be well constrained in the intermediate and soft states. 
However, there are signs for increased A^h towards 0orb = 0 in 
the intermediate state. 

We compared the observed orbital variability of absorption 
with two wind models. A toy model for a focussed line-driven 
wind without clumps ( |Gies & Bolton) 19866) 1 does not describe 
the data well and we attribute the differences to the absence of 
clumping in the model that would lead to the observed strong 
variations of Ah 31 the same orbital phase. Qualitatively match¬ 
ing the observed Ah distributions at a given orbital phase to a 
toy clumpy wind model (|Owocki & Cohen|2006[|Sundqvist et al. 


|2012| l without a focussed wind and with spherical clumps results 
in the best match for models with porosity length, hoa, on the or¬ 
der of the stellar radius, /?,, at orbital phases 0oib ~ 0.5, in agree¬ 
ment with results from (effectively) single O stars ( |Leutenegger| 
|et al.|[20T3| ). Our models exclude both much lower and much 
higher porosities. The discrepancies between the toy clumpy 
wind model and data at (^orb ~ 0 could be due to either a fo¬ 
cussed wind component or a non-spherical shape of the clumps 
or a combination of both. We note here that the existence of a 
focussed wind component in Cyg X-1 is supported by multiple 


lines of evidence such as optical measurements (Gies & Bolton 
|1986bl[So^rs et al.|1998| l. 

While we caution against an over-interpretation of details of 
the presented simple models they firstly demonstate that RXTE 
data of Cyg X-1 can indeed be used for an analysis of the long¬ 
term orbital variability of absorption and secondly clearly de¬ 
mand better models. In particular, a thorough exploration of the 
winds will require a model that includes both, the focussed and 
the clumpy structure of the wind, and possibly also different 
clump shapes and/or sizes. 
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Fig. A.l. Predictions of Nu variability assuming a clumpy wind model with the total number of clumps per flow time N = 1000 and varying 
terminal porosity length /i„ = 10, 3.162, 0.316, and 0.1/?, (see Sect. |4.1.3| for explanation of parameters and Fig.|^for the case of = 1). Left: 
model histograms for the distribution of equivalent absorption column density Nh. Right: average values (circles) and standard deviations (‘error 
bars’ on the average values) in a given orbital phase for the theoretical model (red) and values measured in the hard state of MJD 53900-55375 
(black). Average measured values are shown as ‘average - Ah.ism’ to account for the ISM absorption. 
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Fig. A.2. Predictions of Ah variability assuming a clumpy wind model with the total number of clumps per flow time N and varying terminal 
porosity length /i„. Shown are average values (circles) and standard deviations (‘error bars’ on the average values) in a given orbital phase 
for the theoretical model (red) and values measured in the hard state of MJD 53900-55375 (black). Average measured values are shown as 
‘average - Ah.ism’ to account for the ISM absorption. For the case of A = 1000 see Figs.|^and |A.1| 
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